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ABSTRACT 

The radiative component of the Universe has a characteristic impact on both large 
scale structure (LSS) and the cosmic microwave background radiation (CMB). We 
use the recent WMAP data, together with previous CBI data and 2dF matter power 
spectrum, to constrain the effective number of neutrino species A^e// in a general 
Cosmology. We find that Neff = 4.31 with a 95 per cent C.L. 1.6 < Neff < 7.1. If 
we include the Hq prior from the HST project we find the best fit N^ff — 4.08 and 
1.90 < Neff < 6.62 for 95 per cent C.L. The curvature we derive is still consistent with 
flat, but assuming a flat Universe from the beginning implies a bias toward lower N^f f, 
as well as artificially smaller error bars. Adding the Supernovae constraint doesn't 
improve the result. We analyze and discuss the degeneracies with other parameters, 
and point out that probes of the matter power spectrum on smaller scales, accurate 
independent as measurements, together with better independent measurement of Hq 
would help in breaking the degeneracies. 

Key words: gravitation - cosmology: theory - large-scale structure of Universe - 
cosmic microwave background - dark matter 



1 INTRODUCTION 



The quest for cosmological parameters is the major goal of 
cosmology, and we are experiencing an exciting era in which 
newer and better data allow for the precise determination 
' of an increased number of parameters. In this paper we fo- 
I, cus on one specific parameter, the amount of relativistic en- 
5_j ■ ergy density at recombination/equivalence, and we discuss 
the implications of the current CMB and LSS observations. 
The relativistic energy density is often parametrized by the 
equivalent number of standard model neutrino species N^ff. 
The standard model predicts three neutrino species, plus a 
correction imposed by neutrinos not being completely decou- 
pled during electron-positron annihilation {AN^ff = 0.03) 
(Dolgov et al. 1997, Gnedin & Gnedin 1998) and by the finite 
temperature QED correction to the electromagnetic plasma 
{AN^ff = 0.01) (Lopez et al. 1999, Mangano et al. 2002). 

Departures from the standard model which imply a 
variation in N^ff may be due to decaying dark matter 
(Bonometto & Pierpaoli 1998, Hannestad 1998, Lopez et 
al 1998, Kaplinghat & Turner 2001 or quintessence (Bean 
et al 2001). In these scenarios N^ff can be either smaller or 
greater than 3.04, and the constraints on N^ff implied by 
CMB and LSS need not to agree with the ones implied by 
Big Bang Nucleosyntesis (Kaplinghat & Turner 2001). We 
will analyze and discuss the CMB-I-LSS results on N^f / . Un- 
like previous work on this subject, we allow for a cosmology 
with general curvature. We will use the CMB data from the 



WMAP satellite (Bennet et al. 2003), and complement them 
with the CBI ones (Pearson et al. 2002). As for the LSS, we 
use the 2dF data (Percival et al 2002). We use a general cos- 
mology in our analysis because at the present time we do not 
have an independent probe of the flatness of the Universe 
other than the CMB. It is therefore fictitious and mislead- 
ing to restrict ourselves to flat cosmologies in the so called 
"precision cosmology" era. We will show that allowing for 
non-zero curvature changes quite significantly the results 
on N^ff, despite the fact that the inferred constraints on 
curvature are still compatible with a flat Universe. 

The paper is organized as follows: in section 2 we dis- 
cuss the effects of N^f / on the matter and radiation power 
spectra, in section 3 we discuss the CMB-fLSS constraints, 
while section (4) is dedicated to conclusions. 



2 EFFECTS OF A RELATIVISTIC 

BACKGROUND ON THE MATTER AND 
RADIATION POWER SPECTRA 

An excess of energy in relativistic particles with respect 
to the one predicted by the standard model has an im- 
pact on both the radiation and the matter power spectrum, 
for different reasons. The increased relativistic energy de- 
lays equivalence, and scale entering the horizon at equiva- 
lence is bigger. Since the amplitude of dark matter perturba- 
tions on smaller scales is frozen during the radiation dom- 
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inated era, the matter power spectrum turnover is shifted 
toward lower fc's. The effect can be quantified via the shape 
parameter T ~ nm7i(3*/3.36)^''^ (White, Gelmini & Silk 
1995), where represents the rolativistic degrees of free- 
dom (g* = 2-|-0.454Afe//)- For a given power at large scales, 
the small scale power s therefore reduced, and a lower value 
of f78 is implied. 

An increased radiation density implies a smaller con- 
formal time and therefore a smaller sound horizon. As a 
consequence, the peaks of the CMB power spectrum, which 
reflects the peaks/valleys of the oscillations at recombina- 
tion, arc shifted toward smaller scales (bigger Vs) (Pier- 
paoli & Bnometto 1999). In addition, since recombination 
occurs when the Universe is not completely matter domi- 
nated, the early Integrated Sachs- Wolfe (ISW) effect causes 
an enhanced first peak. 

A description of parameter degeneracies in a flat Uni- 
verse is given by Bowen et al. (2002), where they flnd that 
Nef f is mainly degenerate with flm and the spectral index n. 
They show that almost complete degeneracy can be obtained 
by keeping z^q, Cthh"^ and R (the position of the first acoustic 
peak with respect to a reference model) fixed while varying 
Neff produces almost degenerate power spectra. Here we 
extend the analysis to the real data and consider general 
curvature. 

Allowing for possible curvature affects the radiation 
power spectrum. In particular, the position of the peaks is 
strongly dependent on curvature. A given comoving scale 
at the last scattering surface subtends an angle on the 
sky that is curvature-dependent. If the Universe is open 
(f2fc > 0) [ or closed (Oj; < 0)] then the angle associated 
with a particular scale at last scattering is smaller [larger] 
than the corresponding flat Universe one. Consequently, the 
peaks/valleys of the CMB power spectrum will be located 
at larger [smaller] Vs. A slightly close Universe can compen- 
sate the shift of the peaks toward higher Vs implied by an 
increased radiation background. We therefore expect some 
degeneracy between curvature and N^ff, and if we aim to 
determine both from the CMB measurements, we ought to 
treat both of them as free parameters in the data analysis. 



3 CMB AND LSS CONSTRAINTS 

We used a modified version of the COSMOMC (Lewis & Bri- 
dle 2002) package to compute the likelihoods also including 
the WMAP new results (Verde et al. 2003). We considered 
only CDM adiabatic perturbations, and the following set of 
parameters, with usual definitions: fidm, f^A, ^b, f^fe, Hq, 
Nv, ris, As, Zre (rcionizatiou redshift). In our analysis we 
use the WMAP and CBI data together with the 2dF power 
spectrum, and we discuss the effect of adding other priors. 

First we compare the results for a fiat model (fifc = 0) 
with those for a general Universe when Ho to vary only 
in the range 64: < Hq < 80 corresponding to the la al- 
lowed range of the HST result (Preedman et al. 2001). In 
fig.l we present the likelihoods for each parameter, after 
marginalization over all the others, for fiat and generally 
curved cases. For the number of effective neutrinos we find a 
best fit Neff = 4.08 (iVe// = 3.70) for the general curvature 
(fiat) case, with 2.23 <Neff< 6.13 (1.82 <Neff < 5.74) at 
the 95% C.L. The assumption of a fiat Universe biases the 



result toward low iVe// and shrinks the inferred error bars. 
Our results are in agreement with what found by Crotty et 
al. (2003) with the assumption of fiat geometry. Notice that 
the general curvature analysis is important at the present 
time because of the great improvement that the WMAP re- 
sults have implied on the CMB power spectrum. Constraints 
from CMB and LSS on fiat Universes previous to the WMAP 
were much weaker: N^ff < 15 at 95% C.L. (Hannestad 2001, 
see also Hansen et al., 2002). 

In fig. 1 the fiat and curve case present the similar upper 
limit for N^ff mainly because we imposed a tight top-hat 
prior on Hq (< 80). The likelihood indeed seems to pre- 
fer high values for Ho, which may allow for extra radiative 
component. 

For the general curvature case, we now proceed to dis- 
cuss the effect of various priors on the determination of Neff. 

In fig. 2 we present the results obtained allowing for 
a wider Ho (and Zre) range. The best fit for N^ff in this 
case is Neff = 4.31, with a 95 per cent C.L. range : 1.6 < 
Neff < 7.1 from CMB and 2dF only, and Neff = 4.08 with 
1.9 < Neff < 6.62 when the Hq prior from the HST project 
is included. 

As for the other parameters, we note that the inclusion 
of extra relativistic energy causes a higher Hs values than 
in the standard case (Spergel et al 2003), a higher Zre and 
a slightly closed Universe {Qk = —0.013 ± 0.015). We ar- 
gue that, since a high Neff boosts the first peak through 
the early ISW effect, a higher Us combined with an early 
reionization of the Universe can still ensure a good fit to the 
CMB data. Matter power spectrum data on smaller scales 
than the ones probed by 2dF (e.g. Ly-a forest data) may be 
used in breaking this degeneracy, because the matter power 
spectrum is only sensitive to Us and not to Zre- Data of the 
matter power spectrum at very small scales would be in- 
credibly sensitive to small increases in n^. However, there is 
still much debate on small scale data, since their interpre- 
tation may be complicated by the non-linear growth of the 
fiuctuations. We adopted here a conservative approach and 
chose not to include them in the analysis. In fig. 2 we note 
that the derived value of as is quite high (0.97 ± 0.1) yet 
because of the favoured high Us values. Actual quotations 
of the as value range between 0.75 and 1 ( Pierpaoli et al. 
2003, and references therein), so that a broad prior on as 
would already decrease the allowed Us and improve the er- 
rors on Neff. Consistent 5 per cent accuracy measurements 
of as from both weak lensing and clusters would greatly im- 
prove the constraint on Neff. As an example, in fig. 2 we 
plot the curves obtained with an hypothetical prior on as 
with the typical scaling derived from cluster and weak lens- 
ing: 0-8 (O„/0.3)°-® = 0.85 ± 0.05. Such prior would imply a 
lower Qrn and Neff = 3.8 ± 1.6. 

As for the other parameters, we note that Qm and 
are well constrained in the range: Qm = 0.29 ± 0.06 and 
Oa = 0.72 ± 0.05 (1 a error). 

Hannestad (2003) while analyzing fiat Universes uses a 
tight prior on Qrn derived from the Supernova Cosmology 
Project (Perlmutter e al. 1999). As a consequence he derives 
a small Nef / value with small error bars. We run a separate 
chain treating CMB-|-LSS-|-SN data in a general cosmology 
with COSMOMC, and found that the inclusion of the Su- 
pernovae in the analysis doesn't particularly improve the 
results on Neff (see fig. 2). Including the SN in the anal- 
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Figure 1. The marginalized likelihoods for the parameters un- 
der consideration. We have assumed here a top-hat prior on Hq 
corresponding to the Icr interval allowed by the HST key project 
results. The solid line is for a general curvature, the dotted cor- 
responds to the flat Universe case. The general curvature tends 
to push the constraints on N^ff toward higher values. The same 
upper limits on N^f f are probably due to the upper limit imposed 
on Hq (< 80). Notice that Zre in the range considered here is not 
constrained by the data. 

ysis slightly improves the fi^ and CIa determination, but 
the errors on Neff are practically unchanged. We conclude 
that the strong constraints obtained by Hannestad (2003) 
are mainly due to the restriction to a flat Cosmology and 
to the artificiaUy small prior assumed on f2m as a way of 
introducing the SN constraint. 

In table 1 we plot the correlation matrix for the param- 
eters under consideration. Note that A^^ is mostly correlated 
with Q^rnh^ and flk- A better independent estimate of the 
redshift of equivalence (typically probing flmh), combined 
with a better independent measurement of Hq, may help in 
breaking the degeneracy. 



4 CONCLUSIONS 

The new CMB data (WMAP and CBI) together with the 
matter power spectrum derived by the 2dF galaxy survey 
can constrain the effective number of neutrino species much 
more precisely than previous experiments. Previous esti- 
mates of Neff were derived under the assumption of nul 
curvature. Since we don't have any independent confirma- 
tion of the flatness of the Universe other than the CMB 
itself, we argue that the curvature should be kept as a free 
parameter in the estimation of N^ff. We compare the re- 
sults derived from the two different hypothesis. Applying a 
top hat prior for Hq (64 < Hq < 80), we find for a flat 
Universe 1.82 < N^ff < 5.74 at 95 per cent C.L., with 
a best fit of Neff = 3.70, while with general curvature 
2.23 < Neff < 6.13 with a best fit Neff = 4.08 . Allowing 
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Figure 2. The marginalized likelihoods in the case of a general 
cosmology. The short-dashed line only consider CMB-|-2dF data, 

the dotted includes the Hq prior from the HST project, and the 
solid also includes SN data. Neff is restricted to be < 6.6 at 95 
per cent C.L., and tends to be negative but is still consistent 
with flat. Note the high ris and Zre values. The long— dashed line 
is obtained adding an liyijothctical prior on crs with the typical 
scaling from clusters and weak lensing. 



for general curvature shifts the acceptance range of Neff 
toward higher values mainly because the curvature tends to 
compensate the effect of Neff on the peak locations in the 
CMB power spectrum. 

In the case of general curvature, we explored a wider 
range in Ho and Zre and applied different priors. We find 
1.6 < Neff < 7.1 (best fit Neff = 4.31) at 95% C.L. from 
CMB and 2dF only, and Neff = 4.08 with 1.9 < Neff < 
6.62 when the Hq prior from the HST project is included as 
a proper Gaussian prior. No significant modifications derive 
from the inclusion of the SN constraint. 

By looking at the likelihood distribution for each pa- 
rameter after marginalization over all the others we conclude 
that the inclusion of an extra relativistic component would 
suggest a higher expansion rate Hq , a higher spectral index 
Us and flk slightly negative, if compared to the standard 
analysis with three neutrinos (Spcrgcl ct al. 2003). 

We analyze the correlations between the various param- 
eters and conclude that Neff is most degenerate with Q.dmh'^ 
and fife. We argue that other independent measurements 
of the matter power spectrum, like precise determinations 
of (78 from clusters and lensing or probes at smaller scales 
from the Ly-a forest, would help in constraining the epoch 
of equivalence and therefore would improve the results on 
Neff. Moreover, it would greatly improve the constraints 
on the large Us now allowed. 
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Table 1. Correlation matrix for the parameters in fig 1. We assume here a general curvature, and 64 < Hq < 80. 









Ho 


Zre 






Us 


As 




Age 






1.00 


-0.29 


0.14 


0.39 


0.15 


-0.44 


0.78 


0.70 


0.28 


0.07 


-0.30 


Qw™ h 

^ "am"' 


-0.29 


1.00 


0.25 


-0.07 


-0.10 


0.76 


-0.08 


-0.06 


-0.69 


-0.81 


0.65 


Ho 


0.14 


0.25 


1.00 


0.11 


0.48 


0.17 


0.22 


0.16 


0.46 


-0.75 


-0.57 




0.39 


-0.07 


0.11 


1.00 


-0.38 


0.18 


0.71 


0.87 


0.27 


0.01 


-0.12 




0.15 


-0.10 


0.48 


-0.38 


1.00 


-0.58 


-0.16 


-0.10 


0.17 


-0.28 


-0.47 




-0.44 


0.76 


0.17 


0.18 


-0.58 


1.00 


-0.05 


-0.08 


-0.35 


-0.56 


0.50 


Us 


0.78 


-0.08 


0.22 


0.71 


-0.16 


-0.05 


1.00 


0.89 


0.28 


-0.08 


-0.20 


As 


0.70 


-0.06 


0.16 


0.87 


-0.10 


-0.08 


0.89 


1.00 


0.20 


-0.05 


-0.14 


Ha 


0.28 


-0.69 


0.46 


0.27 


0.17 


-0.35 


0.28 


0.20 


1.00 


0.22 


-0.95 


Age 


0.07 


-0.81 


-0.75 


0.01 


-0.28 


-0.56 


-0.08 


-0.05 


0.22 


1.00 


-0.11 




-0.30 


0.65 


-0.57 


-0.12 


-0.47 


0.50 


-0.20 


-0.14 


-0.95 


-0.11 


1.00 
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